Common quantum phase transition in quasicrystals and heavy-fermion metals 
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Extraordinary new materials named quasicrystals and characterized by noncrystallographic ro- 
tational symmetry and quasiperiodic translational properties have attracted scrutiny. Study of 
quasicrystals may shed light on the most basic notions related to the quantum critical state ob- 
served in heavy-fermion metals. We show that the electronic system of some quasicrystals is located 
at the fermion condensation quantum phase transition without tuning. In that case the quasicrys- 
tals possess the quantum critical state with the non-Fermi liquid behavior which in magnetic fields 
transforms into the Landau Fermi-liquid one. Remarkably, the quantum critical state is robust 
despite the strong disorder experienced by the electrons. We also demonstrate for the first time that 
quasicrystals exhibit the typical scaling behavior of their thermodynamic properties such as the 
magnetic susceptibility, and belong to the famous family of heavy-fermion metals. Our calculated 
thermodynamic properties are in good agreement with recent experimental observations. 
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PACS numbers: 71.23.Ft, 71.27.-f-a, 05.30.Rt 



When encountering exciting behavior of strongly cor- 
related metals, we anticipate to learn more about quan- 
tum critical physics. Such an opportunity may be pro- 
vided by quasicrystals (QCs). These, characterized by 
the absence of translational symmetry in combination 
with good atomic arrangement and rotational symmetry, 
can be viewed as materials located between crystalline 
and disordered solids. QCs, approximants and related 
complex metallic phases reveal very unusual mechanical, 
magnetic, electronic transport and thermodynamic prop- 
erties. The aperiodicity of QCs plays an important role at 
the formation of the properties since the band electronic 
structure governed by the Bloch theorem cannot be well 
defined. As an example, QCs exhibit a high resistivity 
although the density of states (DOS) at the Fermi energy 
is not smalli. One expects transport properties to be de- 
fined by a small diffusivity of electrons which occupy a 
new class of states denoted as "critical states", neither 
being extended nor localized, and making the velocity of 
charge carriers very lowi. Associated with these critical 
states, characterized by an extremely degenerate confined 
wave function, are the so-called "spiky" DOS^i^. These 
predicted DOS are corroborated by experiments reveal- 
ing that single spectra of the local DOS demonstrate a 
spiky DOS^i^. Clearly these spiky states are associated 
with flat bands^ii. On one hand, we expect the proper- 
ties related to the itinerate states governed by the spiky 
states of QCs to coincide with that of heavy-fermion met- 
als, while on the other hand, the pseudo localized states 
may result in those of amorphous materials. Therefore, 
the question of how quasicrystalline order influences the 
electronic properties in quasicrystals, whether these re- 
semble those of heavy-fermion (HF) metals or those of 
amorphous materials, is of crucial importance. 

Recently, experimental measurements on the gold- 



aluminium-ytterbium quasicrystal Au5iAl34Ybi5 with a 
six-dimensional lattice parameter ad = 0.7448 nm have 
revealed a quantum critical behavior with the unusual 
exponent a ~ 0.51 defining the divergency of the mag- 
netic susceptibility x °^ ^ " temperature T 0^. 
The measurements have also exposed that the observed 
non-Fermi liquid (NFL) behavior transforms into Lan- 
dau Fermi liquid (LFL) under the application of a tiny 
magnetic field H , while it exhibits the robustness against 
hydrostatic pressure; the quasicrystal shows also metal- 
lic behavior with the T— dependent part Ap of the re- 
sistivity, Ap cx T, at low temperatures^. All these facts 
challenge theory to explain a quantum criticality of the 
gold-aluminum-ytterbium QC characterized by the un- 
usual exponent and robust against hydrostatic pressure 
but destroyed by tiny magnetic fields. 

In this communication we uncover that a quantum crit- 
ical point of Au5iAl34Ybi5, generating the NFL behav- 
ior, is a fermion condensation quantum phase transition 
(FCQPT)2ii^ and also present the first explanation of the 
low temperature thermodynamics in magnetic fields. We 
explain the robustness of the quantum critical behavior 
against the hydrostatic pressure, and how the applica- 
tion of a weak magnetic field destroys the behavior and 
makes the system transit from the NFL to LFL behav- 
ior. We also demonstrate that there is a general mech- 
anism underlying the NFL behavior of HF metals and 
quasicrystals, leading to a scaling behavior. 

We start with constructing a model to explain the chal- 
lenging behavior of the gold-aluminum-ytterbium QC. 
Taking into account that the spiky states are associated 
to fiat bands^ which are the generic signature of FCQPT, 
we safely assume that the electronic system of the gold- 
aluminum-ytterbium QC Au5iAl34Ybi5 is located very 
near FCQPTS. Thus, Au5iAl34Ybi5 turns out to be lo- 
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cated at FCQPT without tuning this substance with the 
pressure, magnetic field etc. We expect that the system 
exhibits the robustness of its critical behavior against the 
hydrostatic pressure since the hydrostatic pressure does 
not change the topological structure of QC leading to 
the spiky DOS and, correspondingly, flat bands. As we 
will sec, the spiky DOS cannot prevent the field-induced 
Fermi liquid state. 

To study the low temperature thermodynamic and 
scaling behavior, we use the model of homogeneous 
heavy-fermion liquid^. This model avoids the compli- 
cations associated with the anisotropy of solids and con- 
sidering both the thermodynamic properties and NFL 
behavior by calculating the effective mass M*{T,H) as 
a function of temperature T and magnetic field H. To 
study the behavior of the effective mass AI*{T,H), we 
use the Landau equation for the quasiparticle effective 
mass. The only modification is that in our formalism 
the effective mass is no longer constant but depends on 
temperature and magnetic field. For the model of homo- 
geneous HF liquid at finite temperatures and magnetic 
fields, this equation takes the fornt^ i^^i^^ 
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where M is a bare electron mass, F„^ai (pf, p) is the Lan- 
dau interaction, which depends on Fermi momentum pp^ 
momentum p and spin a. Here we use the units where 
h = kB = I- The Landau interaction has the forroi^ 



Fa,a' ip,p') = 



S^E[r 
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where E[n\ is the system energy, which is a functional 
of the quasiparticle distribution function n — liiii^. It can 
be expressed as 



n„{p,T) = <.l + exp 



(e(p,r)-Ata) 



T 



(3) 



where e(p, T) is the single-particle spectrum. In our case, 
the chemical potential fj, depends on the spin due to Zee- 
man splitting = M i ^bH, fJ-B is Bohr magneton. 
The single-particle spectrum is a variational derivative 
of the system energy E[n^{p,T)] with respect to the 
quasiparticle distribution function or occupation num- 
bers n, e{p,T) = SE[n{p)]/6n. In our case F is fixed 
by the condition that the system is situated at FCQPT. 
The variational procedure, being applied to the func- 
tional E[ncr{p,T)], gives using the form for e(p,r) ~ 
£,(p,T) = eK(p,T)], 



dp 



P_ 

M 



E 



Fa,a, (P,Pl) 



5nai(pi,T) (Ppi 



dp 



(4) 



Equations ([3]) and (|4|) constitute the closed set for self- 
consistent determination of ea{p,T) and n„{p^T). The 
sole role of the Landau interaction is to bring the sys- 
tem to FCQPT point, where Af* ^ cx) at T = 
and H ^ 0, and the Fermi surface alters its topol- 
ogy so that the effective mass acquires temperature and 
field dcpcndence^rJiiH. Provided that the Landau in- 
teraction is an analytical function, at the Fermi surface 
the momentum-dependent part of the Landau interac- 
tion can be taken in the form of truncated power series 
F = aq^ + bq^ + cq^ + where q = p — p', a, 6 and c 
arc fitting parameters which are defined by the condition 
that the system is at FCQPT. At this point, the term 
\/M* vanishes and Eq. ([T|) becomes homogeneous. It 
can then be solved analytically^ii^. At = 0, the effec- 
tive mass strongly depends on T demonstrating the NFL 
behavioriSii^ 
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At finite T, the application of magnetic field H drives 
system to the LFL region with 



M*{H) ~ qhH-'^/^. 



(6) 



On the other hand, the interaction F can be repre- 
sented by a nonanalytic function with the first deriva- 
tive 9F(pi, p2)/(?Pi oo at pi ^ p-^^- In that case 
a direct inspection of Eq. ([1]) shows that dX H ~ 0, 
l/M* oc TM*, and we obtain 



M*{T) ~ arT-^'^. 
In the same way, we obtain 

M*{H) ~ ghH-^/^, 



(7) 
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with ax and an are parameters. Taking into account 
the discussed above specific case the gold-aluminium- 
ytterbium QC Au5i Al34Ybi5, we safely assume that the 
Landau interaction is not an analytical function and the 
behavior of M* is given by Eqs. ([7]) and Remem- 
bering that x(2^) 0^ M*{T), we arrive at a good agree- 
ment of the critical exponent defined by Eq. ([7]) with 
the experimental value of a ~ 0.51 ~ 1/2. A few re- 
marks on the transport properties of QC are in order 
here. In calculations of low-temperature resistivity, we 
employ a two-band model, one of which is occupied by 
heavy quasiparticles, with the effective mass given by Eq. 
([7]), while the second band possesses a LFL quasiparti- 
cles with a T— independent effective massif. As a result, 
we find that the quasiparticles width 7 cx T and that 
the T— dependent part of the resistivity Ap cx T. This 
observation is in accordance with experimental fact&S. 

At finite H and T near FCQPT, the solutions of Eq. 
([T]) M*{T,H) can be well approximated by a simple 
universal interpolating function. A deeper insight into 
the behavior of M* (T, H) can be achieved using some 
"internal" scales. Namely, near FCQPT the solutions 
of Eq. ID) exhibit a universal scaling behavior so that 
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M* (T, H) reaches it maximum value M^j at some tem- 
perature Tfnax iJ— li^. It is convenient to introduce the 
internal scales M'^j and T^ax to measure the effective 
mass and temperature, respectively. Thus, we divide the 
effective mass M* and the temperature T by their max- 
imal values, M^^i and Tmax respectively. This generates 
the normalized effective mass — Ad* /M^j and tem- 
perature Tn = T/Tma,^- Near FCQPT the normahzed 
solution of Eq. ([T]) M^{T]sf) with a nonanalytic Landau 
interaction can be well approximated by a simple univer- 
sal interpolating function. The interpolation occurs be- 
tween the LFL {M* (xa + bT^) and NFL (M* oc T'^/^) 
regimes and represents the universal scaling behavior of 
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Here a and b are constants, cq = (1 + C2)/(l -I- ci), ci 
and C2 are fitting parameters. The inset to the left panel 
of Fig. [T] shows the scaling behavior of the normalized 
effective mass. It is seen from the inset, that the common 
width W of the LFL and the transition region W (x T 
vanish as — > since Tmax oc H. In the same way, 
the common width of the NFL and the transition region 
tends to zero as soon as T -> 0. 




field H plays the role of the control parameter, driv- 
ing the system outwards FCQPT that occurs at 7J = 
and T = without tuning since the QC critical state 
is formed by singular density of states^i^i^i^i^. It follows 
from Eq. ([9]) and seen from the left panel of Fig. [U that 
at fixed temperatures the increase of H drives the sys- 
tem along the horizontal arrow from NFL state to LFL 
one. On the contrary, at fixed magnetic field and in- 
creasing temperatures the system transits along the ver- 
tical arrow from LFL state to NFL one. The inset to 
the left panel demonstrates the behavior of the normal- 
ized effective mass versus normalized temperature 
Tn following from Eq. The T~^/^ regime is marked 
as NFL since contrary to the LFL case, where the effec- 
tive mass is constant, the effective mass depends strongly 
on temperature. It is seen that the temperature region 
Ttv ^ 1 signifies a transition regime between the LFL be- 
havior with almost constant effective mass and the NFL 
one, given by T~^/^ dependence. Thus, temperatures 
T ~ Tmax 7 shown by arrows in the inset and the main 
panel, can be regarded as the transition regime between 
LFL and NFL states. The common width W of the LFL 
transition regions ly oc T is shown by the heavy arrow. 
These theoretical results are in good agreement with the 
experimental facts^. The right panel of Fig. [T] illustrates 
the behavior of the dimensionless inverse effective mass 
M/M* versus the dimensionless temperature {T/Tf)^^'^ , 
where Tp is the Fermi temperature of electron gas. To 
calculate M/M*, we use a model Landau functional^ii^ 
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with the nonanalytic Landau interaction 



Magnetic field H 



V{p) = go 



exp(-^o|p| 
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FIG. 1: (color online). Left panel. T — H phase diagram 
of Au5iAl34Ybi5 versus magnetic field H as the control pa- 
rameter. The vertical and horizontal arrows show LFL-NFL 
and NFL-LFL transitions at fixed H and T respectively. At 
H = and T = the system is at FCQPT shown by the 
solid circle. The common width of the LFL and the transi- 
tion regions W oc T are shown by the double arrows. Inset 
shows a schematic plot of the normalized effective mass ver- 
sus the normalized temperature. Transition region, where 
reaches its maximum at T/Tmax = 1, is shown by the 
hatched area. The right panel reports the dimensionless in- 
verse effective mass M/M* versus dimensionless temperature 



{T/Tf 



xl/2 



The line is a linear fit. 



Now we construct the schematic phase diagram of the 
gold-aluminum-ytterbium QC Ausi Al34Ybi5. The phase 
diagram is reported in Fig. [1] left panel. The magnetic 



where the parameters go and /3o are fixed by the require- 
ment that the system be located at FCQPT. It is seen 
from the right panel that the behavior 1/M* oc ^ 
T^/^ is in good agreement with ^ T^-^^ observed 
experimentally^. Such an explanation is consistent with 
the robustness of the exponent 0.51 against the hydro- 
static pressure^ since the robustness is guaranteed by the 
unique singular density of states of QCs^i^i^i^. More- 
over, the crystalline approximant Au5iAl35Ybi4 shows 
the LFL behavior at low temperatures, cx a + bT^-^^ 
with the conventional LFL behavior of the resistivity^. 
We interpret this behavior of through the absence 
of the unique electronic state of QCs, which results in 
the shift of the electronic system of the approximant 
from FCQPT into the LFL region. This shift, creat- 
ing the LFL behavior at T = 0, makes Ti^ax finite 
even at J? = 0. Then, it follows from Eq. ^ that 
1/M* oc x~^ oc a -|- bT^/^ and the approximant is to 
demonstrate the conventional LFL behavior: Ap oc . 
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FIG. 2: (color online). The normalized specific heat {C/T)n 
and normalized magnetic susceptibility xn extracted from 
measurements in magnetic fields H on YbRh2Si2-^ and on 
Au5iAl34Ybig^, respectively. The magnetic fields are given 
in the figure. Our calculations are depicted by the solid curve 
tracing the scaling behavior of (C /T)n = Xn = given by 
Eq. (EJ. 



We now investigate the behavior of x as a function 
of temperature at fixed magnetic fields. The effective 
mass M* (T, H) can be measured in experiments for 
M* (T, H) oc X where x is the ac or dc magnetic suscep- 
tibility. If the corresponding measurements are carried 
out at fixed magnetic field H then, as it follows from Eq. 

X reaches the maximum Xmax at some temperature 
Tmax- Upon normalizing both x a-nd the specific heat 
C/T by their peak values at each field H and the corre- 
sponding temperatures by Tmax, we observe from Eq. Q 
that all the curves merge into a single one, thus demon- 
strating a scaling behavior typical for HE metals^. As 
seen from Eig. [2l xn extracted from measurements on 
Au5iAl34Ybi3^ shows the scaling behavior given by Eq. 
and agrees well with our calculations shown by the 
solid curve over four orders of magnitude in the normal- 
ized temperature. 

In order to validate the phase diagram Eig. [Tl we focus 
on the LEL, NEL and the transition LEL-NEL regions 
exhibited by the QC. To this end, we display in Eig. [3] 
(a) the normalized xn on the double logarithm scale. As 
seen from Fig. [3] (a), xn extracted from the measure- 
ments is not a constant, as would be for a LEL. The two 
regions (the LEL region and NFL one), separated by the 
transition region, as depicted by the hatched area in the 
inset of Fig. [U are clearly seen in Eig. [3] (a) illuminating 
good agreement between the theory and measurements. 
The straight lines in Eig. [3] (a) outline both the LEL 

— 1/2 

and NEL behaviors of xn oc const and xn oc , 
and are in good agreement with the behavior of Af^ dis- 
played in the inset of Eig. [T] In Eig. |31 (b), the solid 
squares denote temperatures T^a.x{H) at which the max- 
ima Xmax of xiT) and, (c), the corresponding values of 
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FIG. 3: (color online), (a): Temperature dependence on the 
double logarithmic scale of the magnetic susceptibility xn 
at different magnetic fields^ shown in the legend. The LFL 
region and NFL one are shown by the solid and dashed arrows, 
respectively. The solid line depicts xn oc T^"'^ behavior, 
(b): The temperatures Tmax at which the maxima of x (see 
Fig. [T]) are located. The solid line represents the function 
Tmax ~ aH, a is a fitting parameter, (c): The maxima Xmax 
versus magnetic field H. The solid curve is approximated by 
Xmax = dH~^^^, see Eq. ([8]), d is a fitting parameter. 



the maxima Xmax(^f) occur. It is seen that the agree- 
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ment between the theory and experiment is good in the 
entire magnetic field domain. It is also seen from Fig. 
|3] (b) that Tmax oc H; thus a tiny magnetic field H de- 
stroys the NFL behavior hereby driving the system to the 
LFL region. This behavior is consistent with the phase 
diagram displaced in Fig. [T] at increasing temperatures 
(T/v — 1) the LFL state first converts into the transition 
one and then disrupts into the NFL state, while at given 
magnetic field H the width W (xT. 

In summary, wc have established for the first time that 
Au5iAl34Ybi5 quasicrystal exhibits the typical scaling 



behavior of its thermodynamic properties, and belongs 
to the famous family of heavy-fermion metals. We have 
also demonstrated that the quantum critical physics of 
the quasicrystal is universal, and emerges regardless of 
the underlying microscopic details of the quasicrystal. 
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